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A  Computer  Program  for  Generating 
Deformation  Mechanism  Maps  for 
Polymers 


1,  Introduction 


The  deformation  mechanism  map  is  a  means  of  representing  the  predominant 
mechanism  of  deformation  within  a  material  over  a  range  of  deformation 
conditions.  M^hanism  maps  have  been  developed  to  describe  creep  behaviour 
in  crystalline  materials  such  a^  metals  and  ceramics  (1, 2],  and  more  recently  to 
describe  deformation  within  polymers  {3-5].  The  program  described  here 
identifies  "ftelds”  within  extension  rate-temperature  space  where  particular 
deformation  mechanisms  operate.  This  program  delineates  these  mechanism 
fields  using  borders,  and  includes  isobars  to  indicate  stress  levels. 

The  program  utilizes  both  theoretical  and  semi-empirical  inodels  of 
deformation  behaviour,  and  while  iir  its  present  form  it  is  restricted  fo 
amorphous  isotropic  polymers,  it  can  be  readily  extended  to  incorporate  other 
polymers,  for  example  semi-crystaliirre,  when  suitable  deformation  models 
become  available. 

The  success  of  a  program  such  as  this  otre  can  be  measured  by  its  ability  to 
simulate  a  purely  empirical  map  constructed  from  experimental  test  data.  The 
deformation  melanism  map  generated  by  this  program  for  PVC  will  be  seen  to 
be  in  very  good  accord  witti  experirtrental  data  in  both  the  predictiori  of 
deformation  mechatusms  and  the  associated  stress  levels.  This  program  has 
been  used  subsequently  as  a  basis  for  gerrerating  rrrechanism  maps  by  Edward 
and  co-workers  [3,  SJ. 


2.  Operation  of  the  Pro^am 


Basically,  the  pro^m  contiains  a  number  of  subroutines  whi(^  comprise 
expressions,  derived  from  thanes  describing  defoirration  within  poijohers  by 
di^erent  meduhisms,  which  relate  the  sta^s  necessary  for  opefatioh  of  each 
ihechanism  to  the  ecte^ioh;  rate  of  ^tehsioir,  temperature  and  oAer  nrateri^ 
propdii«.  At  eadi  value  of  ^rameters  wiUrih  the  extension  rat^ 
tem^rature  s^ce  of  interest;  the  coihputer  cm  rradlly  calculate  the  rnwhanism 
which  requires  the  least  stre^  and  indicate  tMs  on  a  ihap.  A  computer  flow 
chart  des^bing  the  major  features  of  the  program  is  shown  in  Figure  1; 

However,  since  the  moduiiiiof  a  ^lymer  is  ordinarily  dependent  on  both 
time  and  temperature  it  is  hecessaiy  to  Bist  modd  thne  dependencies  so  that 
appropriate  values  ^n  be  used  in  those  equations  which  require  inodulus  data. 
Ihis  program  accepts' a  relaxation  sp^tnim  coiidstihg  of  up  to  5  discfde  sets  of 
relaxation  tittres  and  stfengths.  For.the  pr^ht  exuhple  these  relaxation  times 
and  strengths  are  chosen  empirically  to  ihddei  the  known  time  depeiideiicebf  a 
particular  modulus,  lire  temperature  dependence  then  follows  from  ^suming 
that  the  relaxation  times  obey  the  Williams,  Landel,  Ferry  (WLF)  equation  [6, 7]. 


2,1  Mechanisms  of  Deformation 

A  convenient  temperature  reference  point  for  polymers  is  their  glass  transition 
temperature  since  different  mechanisms  dominate  above  and  below  this 
temperature.  As  there  is  far  less  known  about  the  defontralion  of  polymers 
than  say  metals,  there  are  consequently  orrly  few  theoretical  models  for 
defomration  mechanisms  within  polymers.  This  program  is  therefore  restricted 
to  those  few  fundamental  models  which  are  reasonably  well  established, 
namely;  viscoelasticity,  rubber  elasticity  (T  >  Tg),  Argon  plasticity  (T  <  Tg)  and 
Newtonian  viscosity  (T  >  Tg). 


2,1.1  Rubber  Elasticity,  T>  Tg 

A  model  which  incorporates  the  strain  rate  dependence  of  rubber  eladic 
behaviour  has  been  developed  [8, 9]  through  consideration  of  effective 
crosslinks  (molecular  entanglements)  within  a  polymer  and  the  thermally 
activated  slipping  of  a  molecule  through  such  a  restriction.  The  stress  required 
to  achieve  an  elongation,  X,  at  a  constant  elongation  rate,  X,  is  given  by  [8, 9] 


Og  -  NitT  exp  (Q/itT)(X/X)”‘ 


where  N,  Q  and  m  are  material  constants. 


2.1iZ  Viscoeldsticiiy,  AU  t 

■By  Sfxdfyihg  a  tiiiK,  and  therefore  tem^fature/  dep^dmt  nhc^ulus  of 
ekstidty,  aihechanism  of  d^onhation  to  necessarily  implied.  TOs  being 
\dscoeiastidty.  to  calaUate  the  sinks  wppoftcd  by  tliis.defofmatioh 
meduihism  it  is  hece^aiy  to  substitute  the  modulus  into  ^e  Kltznumh 
superposition  e^tioh  and  evaluate  &is  for  a  renstaht  eiohgatioh  rate  (for 
sniall  extensions). 

Therelevimt  dpiation  is 

t 

o,  -  A.A  JCff-s)* 

O 


where  G(t)  is  determined  through  substitution  of  the  relaxation  spectrum. 


2,13  Argon  Plasticity,  T  <Tg 

This  mechanism,  proposed  by  Argon  [10)>  involves  the  local  rotation  of 
molecular  segments  towards  alignment.with  the  tensile  axis,  this  is  achieved, 
according  to  Argon,  through  th';  introduction  of  kink  pairs  along  the  molecule 
which  enable  it  to  partially  "unfold"  and  «tend  further  in  t)>e  tensile  direction. 
The  stress  required  for  this  mechanism  to  operate  may  be  represented  by  [10] 

-  c,G(1 


2.1,4  Newtonian  Viscosity,  T>Tg 

At  suffidently  warm  temperatures  or  long  times  the  molecules  within  a 
polymer  may  to  some  extent  disentangle  thenaselyes  by  sliding  past  one  another 
in  a  process  termed  "reptation"  [llj.  This  form  of  deforrhation  is  termed 
viscous  flow  and  depends  on  teihperature  through  the  viscosity  term  t)(T). 

Van  Krevelen  [12]  has  shown  that  rdaxation  times  assodated  with  the 
temperature  dependence  of  tj  obey  the  WLF  equation  for  fg<T  <12  Tg,  and 
an  Arrhettius  activation  energy  law  for  T  >  1.2  Tg,  The  stress  can  be  written  in 
terms  of  the  extension  rate  as  follows  112] 


CTn  -  ^A(Ti(r)/3) 


2;2  Input  Data 

The  program  v^tiiies  a  data  ffle  be  crated  in  a  parUoiIar  fonnat  so  that 
Riat^ial  propeitiK  can  be  r&d  m  ^oi%.with  spi^c  data  aardi^hg  to  the 
re^fonm^  of  tte  individi^  medih^ism  When  a  data  value  is  hot 

readied  for  Ae  mechanism  of  interrat,  i£.  m^uliis  reiaxatidn  ^es  for  a 
ciystalltoe  p^l  rahsid^gm  amotphous^  suitable  values 
should  be  substiiuted  in  ^e  in^t  data  in  ttd  appropriate  format  (for  Sample, 
0.0  or  1.0b£:t00). 


2.2^1  GehmlData 


Data  which  are  general  to  the  mechfmism  models  but  specific  to  the  particiiiar 
polymer  Includes  modulus  data  (boA  shear  and  Youngs)  in  both  the  relaxed 
and  unrelaxed  state  (induding  a  spectrum  of  relaxation  times  and  strengths  and 
whether  the  modulus  obeys  an  activation  enei;^  law  or  fhe  WLF  equation),  the 
polymer  melt  and  glass  transition  temperatures,  and  the  percentage  crystallinity 
for  a  semi-crystalline  polymer. 


2.2.2  Mechanism  Specific  Data 

Also  required  for  the  program  is  data  duit  are  essential  to  the  calculations 
performed  for  each  deformation  mechanism.  This  information  indudes  such 
data  as  the  activation  energies  for  spedfic  events,  separation  betw^n  molecular 
kinks,  molecular  radius  and  molecular  frequency  requited  for  Argon  plasUdty 
and  the  glass  transition  temperature  viscosity  for  Newtonian  viscosity. 

Examples  of  input  data  for  PMMA,  PVC,  PET  and  PC  iue  included  in 
Appendix  1. 


2.3  Output  Data  and  Maps 

The  program  is  capable  of  listing  all  the  input  data  in  table  form,  graphing  the 
variation  of  the  moduli  with  both  time  and  deformation  temperature,  creating  a 
table  showing  the  stress  required  by  each  deformation  mechanism  to  operate  at 
each  deformadori  temperature,  and  producing  the  final  mechanism  map 
showing  dondnant  rrtechanisms,  boundaries  between  mkhanisnrs  and  stress 
isobars  in  extension  rate-temperature  space.  An  ^mple  of  computer 
generated  moduli  plots  and  mechanism  maps  for  PMMA  are  ^hown  in 
Appendix  2. 


9 


3.  Correlation  of  Deforniation  Maps  with 
Empiricai  Results 

Ah  eihpiridd'd^orimtion  map  am  be  construe^  ecaminihg  exi^rimental 
data  and  empiri^y  d^ii^g.widdt  dtfpi^tioh  mKhahism  is  most  likely 
opm^  die  ^^t^ion;hxteMion  rate;  stress  wd  tmp^t^ 

values  whidt  this  bq^;  Hpt^  ^  j^ormation  on  a  map  of. extension 
rate  against  tempenituie  (nbimaltebd  with  resj^t  to  the  glass  tr^itioh 
tempeiatute)  generates  im  enipulcd  mediahism  map  whi(^  can  be  compared 
with  that  predicted  by  the  computer  program.  A  compatisoh  of  Figures  2 
snd'3  shows  rearohable  agreement  between  the  empidcal  and  ^eoretic^  maps 
^nd  gives  some  confidence  in  the  iise  of  the  tkebretical  maps  to  ptedid 
deformation  behaviour. 


T/To 


Figure  2:  Empirical  deformation  mechanism  map  showing  experimental  data  points, 
the  logarithmic  stress  leads  associated  with  these  data  points  (e.g.  5  =  l(f  Pa,  with 
solid  litia  representing  boundaries  for  order  of  magnitude  stress  increases)  and  the 
border  betxeeen  mechanism  fields  (*).  Reproduced  from  reference  13]. 
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ligure  3:  Theoretical  deformation  mechanism  map  showing  logarithmic  stress  levels 
(e.g.  4  =  10*  Pa,  and  solid  lines  represent  boundaries  for  order  of  magnitude  stress 
«  increases)  and  the  border  between  mechanism  fields  M.  Reproduced  from 

I  reference  [3], 


j  4.  Conclusion 

i 

I  A  program  has  been  presented  that  Is  based  on  various  theoretical  models  for 

polymer  deformation  and  is  able  to  dedtire  which  of  the  ihechanisnas  proposed 
by  the  deformation  models  is  most  likely  to  operate  at  any  given  combination 
of  extension  rate  and  temperature.  The  program  then  generates  a  deformation 
mechanism  map  in  extension  rate  temperature  space  indicating  the  fields  in 
which  each  mechanism  dorhinates  and  the  str^s  levels  at  which  these 
mechanisms  operate. 

The  deformation  mechanism  map  so  produced  for  PVC  is  seen  to  closely 
resetirble  a  similar  map  deduced  soldy  fiom  experimental  observatiotrs  giving 
some  coirfidence  in  >he  predictive  capabilities  of  the  program. 
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7.  Nomenclature 

elongation 
elongation  rate 
temperature 
Boltzmann's  constant 
number 

activation  energy  for  entanglement  slipping 

strain  rate  sensitivity 

shear  modulus 

time 

viscosity 

constants 


Appendix  1 


Con^uter  input  data  for  a  number  of  conunon  polymers  sourced  from 
experiment  ifod  the  quoted  references.  The  tdininologjr  used  is  fully  described 
in  the  conq>uter  program. 


Polymethylmcthaciylate  (FMMA) 


SHMOD  U,R,  e  NLNKSR 

1.10E>09 

S.OOE-'OS 

1.25Et21 

YMOD  U,R  t  XLNXY 

3.20E«09 

l.SOEtOe 

o.oosm 

SHMOD  UC,RC  «  PHI 

4.30Et09 

1.50E^09 

0.30E+00 

YMOD  UC.RC  C  PSI 

9.00E+09 

4.50Ei09 

O.EOEtOO 

GHEI.HECRIT.KHI 

5.20E+02 

9.00E' 

-02 

O.SOE^-OO 

LOG (RELXN . TIMES) AMOR 

-1.00 

0.00 

1.00 

0.00 

0.00 

LOG (RELXN . TIMES) CRYS 

0.00 

0.00 

0.00 

-3.87 

-2.67 

REfAXN.STRGTH  (AMOR) 

O.ZSE-fOO 

O.SOEtOO 

0.2SE400 

O.OOE^OO 

O.OOE-fOO 

RELAXN.STRGTH (CRYS) 

0.00B4^00 

O.OOE^OO 

O.OOE+00 

0.60E^00 

3.406400 

ACT. ENERGIES (AMOR) 

O.OOE-09 

O.OOE-20 

O.OOE-OO 

O.OOE-OO 

O.OOB-09 

ACT. ENERGIES (CRYS) 

O.OOE-00 

O.OOE-OO 

O.OOE-OO 

1.25E-19 

1.2SE-19 

OBEYS  Q  OR  HLP(AM) 

1.0 

1.0 

1.0 

1.0 

1.0 

OBEYS  0  OR  WLF(CR) 

2.0 

2.0 

2.0 

2.0 

2.0 

TMELT (K) , TGLASS (K) TREF 

4.50Et02 

3.78E+02 

3.93E+02 

EACTl.EMMl.ENNl 

S.60E-2D 

1.30E' 

-01 

O.00E400 

KNKSEP, MOLRAD, HOLPBN 

2.85E-09 

2.8SE' 

-10 

1.00E400 

OMEGA, MOLFRQ.VOLDHS 

SiSOE-fOO 

1.80E+11 

2.00E428 

VaCT1,KKSRT,KINTSI 

2.0E-16 

1 

5.0 

1.20E-23 

ETATG,ETAACT,GHE2 

l.OOE+03 

5.00E-21 

O.OOE-OO 

Plasticized  Polyvinylchloride  (PVC) 


SHMOD  U,R,  8  NlNKSR 

1.00Ef09 

7.00E+05 

1.25E+21 

YMOD  U,R  8  XLNTY 

3.00E+09 

7.00E+05 

3.00E+01 

SHMOD  UC,RC  4  PHI 

3.30E409 

7.00E+07 

3.00E-01 

YMOD  UC,RC  (  PSI 

9.90E+09 

7.00E+07 

3.00E-01 

GHB1,GHE2,KHI 

l.OOE-01 

2.00E-02 

5.00E-01 

LOG  (RELXN .  TIMES)  AMOR 

-5.0 

-2.0 

-1.0 

0.0 

1.0 

LOG  (REUN .  TIMES)  CRYS 

-2.0 

-3.0 

-4.0 

5.0 

6.0 

RELAXN . STRGTK (AMOR) 

O.SOBtOO 

0.50E+00 

O.O0E4OO 

O.O0E4OO 

O.OOE+00 

RELAXN.STRGTH(CRYS) 

0.30P.400 

0.20E+00 

0.50B+00 

o.eoE400 

O.00E400 

ACT. ENERGIES (AMOR) 

1.50E-19 

1.60E-20 

l.OOE-10 

O.OOE-OO 

1.50E-19 

ACT. ENERGIES (CRYS) 

l.OOE-01 

3.10E-18 

3.20E-19 

1.25E-19 

3.30E-20 

OBEYS  Q  OR  WLF(AM) 

1.0 

1.0 

1.0 

1.0 

1.0 

OBEYS  Q  OR  WLF(CR) 

2.0 

1.0 

1.0 

2.0 

2.0 

TMELT (K) , TGLASS (K) TREF 

S.O0E402 

2.70E+02 

3.00E402 

EACTl.EHMl.ENNl 

5.60E-20 

1.30E-01 

O.OOE400 

KNKSEP, MOLRAD, MOLPBN 

2.85E-09 

2.85E-10 

l.OOEtOO 

OMEGA, MOLFRQ, VOLDNS 

5.60E400 

1.80E+11 

2.00E428 

VACT1,KKSRT,KINIST 

2.00E-16 

e.OOEiOO 

1.20B-23 
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Folyeftylehe  teiephthalate  (FCT) 


SHHOD  UiR,  e  NLNKSR 
KMOD  U,R  8  XLNTY 

1.10E409 

3.00B'f09 

7.00E+05 

T.dOE-fOS 

1.25E+21 

3.00E+01 

LOG  (^U(N ;  flHES)  AHOR 

-3id 

-2.0 

-1.0 

0.0 

1.0 

LOG  (RELXNiTIffiS)CRYS 

2'.0 

3.0 

4.0 

5.0 

6.0 

rslaxn;strgih(ah6r) 

O.SOEiOO 

O.SOE+OO 

d.OOE+OO 

O.dOE-fOO 

O.OOE^OO 

RELAXN.SIRGIH(CRYS) 

O.OOE-t-dO 

b.odB^-oo 

O.dOE+OO 

OiOOE-fOO 

O.OOE^-OO 

ACT. ENERGIES (AHOR) 

1.50E-19 

l.S0Er20 

l.OOE-10 

I.SOE-IS 

1.50E-19 

ACT. ENERGIES (CRYS) 

l.OOE-01 

3.I0E-18 

3.20E-19 

3.00B-18 

3.30E-20 

OBEYS  Q  OR  WLF(M) 

i.o 

1.0 

1.0 

1.0 

1.0 

OBEYS  Q  OR  WLF(CR) 

2.0 

1.0 

1.0 

2.0 

2.0 

TMELT  (K) ,  TGLASS  (K)  FREQG 
EACT1,EHH1,ENN1 

KNKSEP , MOLRAD , HOLPBN 

ONEGA, HOLFRQ.VOLONS 
VACTl.KKSRT.KINTST 

5.45E+02 

5.60E-20 

2.8SE-09 

S.eOE-fOO 

2.00E-16 

3.45Et02 

1.3dE-01 

2.85E-10 

1.80E+11 

6.0 

l.OOE+00 

O.OOE-fOO 

l.OOE+dO 

2.00E+28 

1.20E-23 

Polycarbonate  (PC) 


SHMOD  U,R,  8  NLNKSR 
YMOD  U,R  8  XLNTY 

l.lOB+09 

3.12Et09 

7.00E+05 

7.00E+05 

1.2SE-f21 

3.00E+01 

LOG (RELXN .TIMES) AMOR 

-3.0 

-2.0 

-1.0 

0.0 

1.0 

LOG(RELXN.TIMBS)CRYS 

2.0 

3.0 

4.0 

S.O 

6.0 

i, 

RELAXN . STRGTH (AMOR) 

O.SOE+OO 

O.SOEtOO 

0.00E+C7 

O.OOE-fOO 

O.OOE-fOO 

RELAXN.SIRGTH(CRYS) 

O.OOEiOO 

O.OOE-fOO 

O.OOE'^OO 

O.OOE-fOO 

O.OOE-fOO 

i 

\ 
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Appendix  2 

i 

Computer  generated  plots  of  shear  modulus  (G)  and  Poissons  ratio  (P)  versus 
(i)  normalized  temperature  (T/Tg)  and  (ii)  time;  and  computer  generated 
mechanism  maps  showing  the  operation  of  viscoelastiaty  (V),  rubber 
elasticity  (E),  Argon  plasticity  (A)  and  Newtonian  viscosity  (N)  within  PMMA 
over  extension  rate  -  temperature  for  increasing  values  of  extension  (X). 

Asterisks  (*)  represent  borders  between  mechanisms  and  CTOSses  (X) 
interspersed  with  numerals  represent  order  of  magnitude  increases  in  stress 
magnitude,  i.e.  3  =  10^  Pa. 
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